Nanodroplets on rough hydrophilic and hydrophobic surfaces 
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We present results of Molecular Dynamics (MD) calculations on the behavior of liquid nan- 
odroplets on rough hydrophobic and hydrophilic solid surfaces. On hydrophobic surfaces, the con- 
tact angle for nanodroplets depends strongly on the root mean square roughness amplitude, but it 
is nearly independent of the fractal dimension of the surface. Since increasing the fractal dimension 
increases the short- wavelength roughness, while the long- wavelength roughness is almost unchanged, 
we conclude that for hydrophobic interactions the short-wavelength (atomistic) roughness is not very 
important. We show that the nanodroplet is in a Cassie-like state. For rough hydrophobic surfaces, 
there is no contact angle hysteresis due to strong thermal fluctuations, which occur at the liquid- 
solid interface on the nanoscale. On hydrophilic surfaces, however, there is strong contact angle 
hysteresis due to higher energy barrier. These findings may be very important for the development 
of artificially biomimetic superhydrophobic surfaces. 
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I. INTRODUCTION 

In the year of 1805, Thomas Young and Pierre Si- 
mon de Laplace proposed that an interface between two 
materials has specific energy, the so-called interfacial 
energy, which is proportional to the interfacial surface 
areajll H HI- This concept is the basis for the field of 
wetting, which has become an extremely hot topic in the 
last two decades 0,(1], thanks to biological and high-tech 
applications, ranging from self-cleaning surfaces, micro- 
electronics and thin film coatings, to image formation 
that involve the spreading of liquids on solid surfaces. 

Wetting describes the contact between a fluid and a 
solid surface. Liquids with high surface tension (usu- 
ally reflecting strong intra- molecular bonds), or liquids 
on low-energy solid surfaces, usually form nearly (com- 
plete) spherical droplets, whereas liquids with low surface 
tension, or liquids on high-energy surfaces, usually spread 
out on (or wet) the surfaces. This phenomenon is a result 
of minimization of interfacial energy. Thus, if a surface 
has a high free energy, most liquids will spread on the 
surface since this will usually lower the free energy. 

Wetting phenomena have been widely studied both 
theoretically^, 0] and experimentally]!, Q in connection 
with the physics of surfaces and interfaces. The behavior 
of liquids on smooth solid surfaces is rather well under- 
stood. However, for rough solid surfaces the situation is 
much less clear, even though roughness occurs on prac- 
tically all real surfaces of engineering or biological inter- 
est. Studies (and classification) of disordered and inho- 
mogeneous surfaces [l(| should have significant impact on 
the problem of li quid contact angle and wetting of rough 
substrates [H El 11 EI El 

The fascinating water repellents of many biological sur- 
faces, in particular plant leaves, have recently attracted 
great interest for fundamental research as well as practi- 
cal applications [H HE El El M, M, HI H3 . The abil- 
ity of these surfaces to eliminate water beads completely 




FIG. 1: A droplet on a superhydrophobic surface: The 
droplet touch the leaf only at a few points and forms a ball. 
It completely rolls off at the slightest declination. Adapted 
from Ref. with permission. 



and thereby wash off contamination very effectively has 
been termed the Lotus effect, although it is observed not 
only on the leaves of the Lotus plant (Fig. [1]), but also 
on many other plants such as strawberry, raspberry and 
so on. Water repellents are very important in many in- 
dustrial and biological processes, such as prevention of 
the adhesion of snow, rain drops and fog to antennas, 
self-cleaning windows and traffic indicators, low-friction 
surfaces and cell mobility 0, [25l. [26j. 

Most leaves that exhibit strong hydrophobicity have 
hierarchical surface roughness with micro- and nanos- 
tructures made of unwettable wax crystals, which maxi- 
mize the contact angle with water and most other liquids. 
Fig. [5] shows epidermal cells (microscale roughness) cov- 
ered with wax crystals (nanoscale roughness). The wax 
crystals exhibit a relative high contact angle with wa- 
ter, which is enhanced by the surface roughness. Water 
droplets on the rough wax surface tend to minimize the 
contact between the surface and the droplet by forming 
nearly spherical droplets, as approximately described by 
the two classical models due to Wenzel|27l and Cassiel2a 
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FIG. 2: A leaf surface with roughness on several length scales 
optimized via nature selection for hydrophobicity and self- 
cleaning. Through the combination of microstructure (cells) 
and nanostructure (wax crystals) the macroscopic water con- 
tact angle do is maximized. Adapted from Ref. pij with per- 
mission. 

(see below). As a result the leaves have also a self- 
cleaning property: because of the small adhesion energy 
(and small contact area) between contamination parti- 
cles and the rough leaf 17], during raining water drops 
roll away removing the contamination particles from the 
leaf surface. 

The hydrophobicity of solid surfaces is determined 
by both the chemical composition and the geometrical 
micro- or nanostructure of the surface [{j, l29l. [30j . Under- 
standing the wetting of corrugated and porous surfaces is 
a problem of l ong standing interest in areas ranging from 
textile science [31j to catalytic reaction engineering [32| . 
Renewed interest in this problem has been generated by 
the discoveries of surfaces with small scale corrugations 
that exhibit very large contact angles for water and other 
liquids — in some cases the contact angle is close to 180°. 
Such surfaces are referred to as superhydrophobic[33]. 

In this paper we present results of Molecular Dynamics 
(MD) calculations on the behavior of liquid nanodroplets 
on rough hydrophilic and hydrophobic solid surfaces. We 
find that for hydrophobic surfaces, the contact angle for 
nanodroplets depends strongly on the root mean square 
surface roughness amplitude, but is nearly independent of 
the fractal dimension Df of the surface. For hydropho- 
bic rough surfaces we do not detect any contact angle 
hysteresis. Both results can be explained by the strong 
thermal fluctuations which occur at the liquid-solid inter- 
face on the nanoscale. On hydrophilic surfaces, however, 
strong contact angle hysteresis has been found due to the 
higher energy barrier for interfacial liquid density fluctua- 
tions. These findings may be crucial for the development 
of artificial biomimetic superhydrophobic surfaces. 

II. THEORETICAL BACKGROUND 

In this section we briefly describe some results from 
the theory of the liquid-solid contact angle, which are 
necessary for the interpretation of the numerical results 
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FIG. 3: Liquid droplet on flat substrate. The contact angle 
9 is between (complete wetting) and it. 

presented in Sec. IIVI We emphasize the importance 
of thermal fluctuations for the contact dynamics at the 
nanoscale as compared to micrometer or macroscopic di- 
mensions. 



A. Flat surfaces 

If gravitational effects can be neglected, a liquid 
droplet on a flat substrate forms a spherical cap, see 
Fig. [31 The contact angle 8 is determined by the mini- 
mization of the free energy and depends on the interfacial 
free energies per unit area: solid/liquid 7 s i, solid/vapor 
7sv and liquid/vapor 7i v . Minimizing of the surface free 
energy, with the constrain of fixed volume of the droplet, 
gives the Young's equation, first proposed by Thomas 
Young about two hundred years ago: 

7sl + 7i v cos6» = 7 SV (1) 

Complete wetting corresponds to 9 — 0, and typically 
happens for liquids with low surface tension j\ v , and on 
solids with high surface energy 7 SV . Liquids with high 
surface tension on surfaces with low surface energy tend 
to form droplets with high contact angle 9. Eq. ([T]) 
was deduced for a substrate which is assumed to be per- 
fectly smooth, homogeneous, and rigid. However, in re- 
ality, structured or rough surfaces are quite common. So 
it's necessary to know how the contact angle behaves on 
rough surfaces. 



B. Rough surfaces: minimum free energy state 

Most surfaces of practical interest have roughness on 
many different length scales. For simple periodic surface 
profiles one may develop accurate analytical treatments 
of the liquid droplet contact angle (see e.g., Ref. 
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FIG. 4: Liquid droplet on a rough substrate. At the lowest 
magnification £ the surface appears flat and the liquid contact 
angle is 6q. At increasing magnification surface roughness is 
observed and the liquid will in general only make contact with 
the substrate in some asperity contact regions. 

but for randomly rough surfaces the situation is much 
more complex. For surfaces with random roughness, e.g., 
sclf-affine fractal surfaces (see below), one may develop 
a general theory based on the study of the system at 
different magnifications £, see Fig. @] Here ( = D/X 
where D is the diameter of the droplet-substrate (appar- 
ent) contact area and A the resolution. One can intro- 
duce effective interfacial liquid-solid and solid-vapor free 
energies (per unit area) 7 s i(C) and 7 SV (C) which depend 
on the magnification [351] . At the highest magnification 
£i , corresponding to nanometer (or atomistic) resolution, 
these quantities reduce to those for the flat surface [36|, 
7si(Ci) — 7si an d 7sv(Ci) — 7sv- Since the substrate ap- 
pears flat at the lowest magnification £ = 1, the macro- 
scopic contact angle (corresponding to £ = 1) is obtained 
using the Young's equation with 7 s i and 7 SV replaced by 
7si(l) and 7sv (l), i.e. 

7 s i(l) +7i v cos6»o =7sv(l) (2) 

The change in the surface free energy (per unit area) 
when a liquid with a flat surface is brought in contact 
with the substrate is 

AF/A Q = 7sl (l) - 7sv (l) - 7i v = - 7lv (l + cos6» ) 

where Aq is the (projected) surface area. Note that in- 
creasing contact angle 9q corresponds to a increasing in- 
terfacial free energy. Thus, if a liquid drop can occur in 
several metastable states on a surface, the state with the 
smallest contact angle corresponds to the (stable) mini- 
mal free-energy state. 

Using Eq. © it is trivial to derive the results of the so 
called Wenzei|27j and Cassie 28] models. In the Wenzel 



model it is assumed that complete contact occurs at the 
liquid-solid interface. Thus 

7sv(i) = nfev(Ci) . 7si(i) = nfei(Ci) , (3) 

where r = A/Aq > is the ratio between the surface area 
A of the rough substrate, and the projected (or nominal) 
surface area Aq. Substituting ^ into © gives the con- 
tact angle 9q on the rough surface in terms of the contact 
angle 9 on the microscopically flat surface of the same 
material (Wenzel equation): 

cos 9o — r cos 9 . (4) 

In the Cassie model[28| it is assumed that some air 
(or vapor) remains trapped between the drop and the 
cavities of the rough surface. In this case the interface 
free energy 

7sv(l) = r 7s v(C0, (5) 

7sl (l) = 0r 7sl (Ci) + (1 - 0)(r 7sv (Ci) +7iv), (6) 

where 4> is the fraction of the (projected) area where the 
liquid is in contact with the solid. Substituting ([5]) and 
© in © gives 

cos O = rcos6> - (1 - <p)(l + rcos9) . (7) 

Note that for </> = 1, ([7]) reduces to ((4]). In the original 
Cassie model it was assumed that r = 1. We note that 
while the Wenzel theory is exact if the liquid is in con- 
tact with the substrate everywhere within the nominal 
liquid-substrate contact area, the Cassie theory is always 
approximate and often not very accurate. This is easily 
understood from Fig. [5] which shows the interface be- 
tween a liquid and a solid. </> < 1 is the ratio between the 
projected liquid-solid contact area and the nominal (or 
apparent) contact area A . Because the solid surface is 
curved, the actual liquid-solid contact area will be A^tps 
where s > 1. Analogously, since in general the liquid- 
vapor interface is curved (in spite of the fact that the 
total curvature l/i?i + 1 /R2 may vanish) and tilted (rel- 
ative to the average surface plane), the total liquid- vapor 
interface area is Aq(1 — 4>)s', with s' > 1. Similarly, the 
solid- vapor interface area equals Ao(l — <fi)s" with s" > 1. 
In deriving ([7]) it is assumed that s — s" — r and s' = 1 . 

Of the two states, Cassie and Wenzel, the stable one, 
that is the one with lower free energy, is the one with 
larger cos 9 - Comparing ([4]) and ([7]) shows that the value 
of cos #0 for the Cassie state is larger if 1 + r cos 9 < or 

cos0<-l/r. (8) 

Since r is a measure of the magnitude of the surface 
roughness, we may qualitatively state that only for hy- 
drophobic surfaces (with cos6> < or 9 > 90°) with large 
enough roughness (i.e., large enough r = A/Aq) will the 
Cassie state be the thermodynamically stable state. 
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FIG. 5: The interface between liquid and solid. <f> < 1 is 
the ratio between the projected liquid-solid contact area and 
the nominal (or apparent) contact area Aq. Because the solid 
surface is curved, the actual liquid-solid contact area will be 
Ao(j>s where s > 1. Similarly, since in general the liquid- 
vapor interface is curved (in spite of the fact that the total 
curvature 1/Ri + 1 /R2 may vanish) and tilted (relative to the 
average surface plane) the total liquid-vapor interface area is 
Aq(1 — (f>)s , with s > 1. Similarly, the solid-vapor interface 
area equals Ao(l — 4>)s" with s" > 1. 

The approach described above, where the interface is 
studied at different magnifications, is very general and 
a similar approach has recently been developed for the 
contact mechanics between elastic solids with randomly 
rough surfaces [13] (see also Ref. [HI]). 

It is well known that the roughness of a hydrophobic 
solid (with 9 > 90° on the flat substrate) enhances its 
hydrophobicity. If the contact angle of water on such 
flat solids is of the order of 100° to 120°, on a rough 
or microtextured surface it may be as high as 150° to 
175° [II, 11 111 . Both the Wenzel model and the Cassie 
model can explain this effect. 

Let us consider the simplest surface roughness consist- 
ing of a periodic rectangular roughness profile as illus- 
trated in Fig. [6] (a) (xz- plane) . The free energy (per unit 
surface area) for the Cassie state shown in the figure is 

7C = [(a + 2h)j sv + a-yu + bj s i] /(a + b) 

The free energy for the Wenzel state (complete contact) 
is 

7w = (a + b + 2h)j sl /{a + b) 

Using ([1]) we can write the difference in free energy 

7c - 7w = 7iv[a(l + cos 6) + 2ft,cos 0]/(a + b) 

Thus, the Cassie state has a lower free energy than the 
Wenzel state if 

cos9<-(l + — ) (9) 



FIG. 6: (a) Liquid drop (in the Cassie state) in contact with 
a surface with periodic surface roughness, (b) Even if the 
Cassie state (incomplete liquid-solid contact) is the ground 
state, with an applied pressure p one can squeeze the droplet 
into the Wenzel state. 

which is satisfied only if for the flat surface 6 > 90°, and 
if the ratio h/a is large enough. Note that in this case 
r = A/Aq = (a + b + 2h)/(a + b) so the (approximate) 
criteria (O reduces to 



cos 6 < — 




which is of similar general form as © ■ In Nature strongly 
hydrophobic surfaces are often obtained by covering the 
surface with thin, long (so that h/a ^> 1) hydrophobic 
fibers. Thus, insects which move on top of water, e.g., 
water spiders (see Fig. [7]) have a high density of thin wax 
coated hair on their legs. In addition, the hair fibers have 
nanoscale roughness which traps air and enhances the 
hvdrophobicitv [io] . In this case the water-leg contact will 
be in the Cassie state even when the insect is squeezed 
towards the water by the weight of the insect. 

C. Rough surfaces: activation barriers and 
hysteresis 

Consider a cylindrical cavity as in Fig. [S] (b) and as- 
sume first the Cassie state as in the figure. Let us apply 
a pressure p to the droplet. In this case the liquid will 



FIG. 7: Water spiders have thin hydrophobic (wax coated) 
hair with nanoscale roughness which trap air and enhance 
hydrophobicity. 

bend inwards in the cavity and if the applied pressure is 
larger than a critical value p c , the liquid will be squeezed 
into the cavity (we assume that the air in the cavity can 
leave the cavity, e.g., diffuse into the liquid). It is easy 
to show that the pressure 

p c = -2 llv cos9/R (10) 

where R is the radius of the cavity. To prove this relation, 
note that the pressure work to squeeze the liquid a dis- 
tance h! into the cavity (see Fig. [5]) is given by p c irR 2 h' 
and this must equal the change in interfacial free energy 
which equals 2tt Rh' (j s \ — j sv ) . Using these equations and 
© gives CEU]). 

From HUl) it follows that if 9 < 90° (hydrophilic in- 
teraction), p c < and the liquid will be spontaneously 
sucked into the cavity and will fill out the cavity. If 
9 > 90° (hydrophobic interaction), p c > and for 
nanometer sized cavities, the pressure p c ~ 100 MPa, so 
very high pressures are necessary for squeezing the liquid 
into narrow cavities. However, if the liquid is squeezed 
into the cavity and completely fills the cavity, then the 
resulting Wenzel state is (at least) metastable. How- 
ever, for nanometer sized cavities thermal fluctuations 
may give rise to strong local fluctuations between the 
Cassie (empty cavity) and Wenzel (filled cavity) states. 
This is easy to understand since the energetic barrier 
(for a hydrophobic system) for going from the Cassie 
state to the Wenzel state will be of order e ~ p c iiR 2 h = 
— 2ir RhjivCOsO , and strong fluctuations on macroscopic 
time scales will occur as long as e ~ 0.7 eV or less, and 
strong fluctuations on the nanosecond time scale occur 
if e ~ 0.4 eV or less [note: the rate to jump over a 
barrier of height e is w = v exp(— e/k^T) where typi- 
cally v as 10 12 s _1 ; at room temperature it « 1 s _1 if 
e k, 0.7 eV and w « 10 9 s^ 1 if e w 0.4 eV]. In a typ- 
ical case this condition is satisfied as long as R and h 
are of order of one nanometer or less. In our computer 



FIG. 8: The fraction P of the surface area where the abso- 
lute value of the slope is smaller than |tan#| as a function of 
|tan#|/£. See text for details. 



simulations we do indeed observe very strong thermal 
fluctuations at the liquid-solid interface, in particular for 
rough hydrophobic surfaces, see Sec. IIVB II 

The Wenzel droplets are highly pinned, and the tran- 
sition from the Cassie to the Wenzel state results in the 
loss of the anti-adhesive properties generally associated 
with super hydrophobicity. However, for nanodroplets 
on rough hydrophobic surfaces, we find that the Wenzel 
state is unstable: if the droplet is pressed into complete 
contact with the substrate (Wenzel-like state) and then 
let free, it quickly jumps back to the Cassie-like state due 
to strong thermal fluctuations. For a macroscopic droplet 
on surfaces with long wavelength roughness, the energetic 
barrier towards flipping from the Wenzel to the Cassie 
state may be so large that, even if the Cassie state is 
the minimum free energy configuration, the system may 
remains trapped in the (metastable) Wenzel state for all 
time periods of physical relevance. 



D. Cassie and Wenzel states for randomly rough 
surfaces 

In this section we discuss under which condition one 
expects the Cassie state or the Wenzel state to prevail. 
Consider a rough surface and let z = h(x) be the height 
of the surface at the point x = (x, y). A randomly rough 
surface can be obtained by adding plane waves with ran- 
dom phases: 

/i(x) = J2 S(q)e l[q - x+0(q)1 
q 

where <j>(q) are independent random variables, uniformly 
distributed in the interval [0, 2tt[, and with B(q) = 
(27r/L)[C(q)] 1 / 2 , where L = Al' 2 is the linear size of 



the surface. The surface roughness power spectrum: 



c(q) = ^hf\ Sx <M*)MO)>e <q ' x (11) 

Here h{x) is the surface height profile and (■ ■ •) stands 
for ensemble average. We have assumed that (h(x)) = 0. 
We assume that the statistical properties of the rough 
surface are isotropic, so that C(q) only depends on the 
magnitude q = |q| of the wave vector q. 

For randomly rough surfaces the normalized surface 
area r = A/Aq is given by (see Appendix): 

r= dx (l + xf) 1/2 e- x (12) 
Jo 

where 

e = J d 2 q q 2 C(q) = 2ir jH dq q 3 C(q) = ((\7h) 2 ) (13) 

is the square of the average slope. 

The fraction of the surface where the surface slope s < 
so is given by (see Appendix): 

P(s ) - l-e-( s °/«) 2 . 

Note that as £ — > (corresponding to a flat surface) 
P(sq) — > 1 which is expected because the slope of a flat 
surface is zero and hence smaller than any finite value 
so- Assume that a liquid exhibits the contact angle 9 on 
the perfectly flat substrate. The fraction of the surface 
where the slope |V/i(x)| < |tan0| is given by 

P(tan0) = 1 -e" (tane/ « )2 . (14) 

This function is shown in Fig. [5] If we assume that the 
liquid surface in the liquid-solid non-contact region is flat 
and parallel to the average substrate surface plane, then 
we expect the liquid to only occupy the region where 
the slope is smaller than tan9. Note that more than 
90% of the surface area will have a slope below |tan#| if 
tan#|/£ > 1.5 and in this case the Wenzel state will tend 
to prevail, while more than 90% of the surface will have 
a slope above |tan0| if |tan#|/£ < 0.3 and in this case the 
Cassie state will tend to prevail. For the system we study 
below £ < 2 (see Fig. [T2")) and for the hydrophobic system 
9 w 103° we get |tan0|/£ > 2.2. Thus, one would expect 
the Wenzel state to prevail. However, the numerical data 
(see below) tend to suggest that the system is in a Cassie- 
like state. We attribute this to the strong fluctuations at 
the liquid-solid contact which occurs at the nanoscale, 
which are particularly important for nanoscale droplets. 

III. SIMULATION METHOD 

We have used Molecular Dynamics (MD) to study the 
contact angle and contact angle hysteresis. Here we 
briefly describe the system we studied and how we gen- 
erated the rough substrate surfaces. 




FIG. 9: 3D side-view snapshot of a octane liquid droplet 
on a hydrophobic and rough substrate. The rigid substrate 
comprises 200 x 30 atoms disposed on a square lattice with 
lattice constant a = 2.53 A. These atoms have been randomly 
displaced along the z coordinate, orthogonal to the wall, so 
to reproduce the desired roughness. The Lennard- Jones solid- 
liquid interaction potential V(r) = 4e[(ro/r) 12 — (ro/r) 6 ] with 
ro = 3.28 A, e = 4 meV for hydrophobic substrate and e = 
8 meV for hydrophilic substrate. 



A. Molecular dynamics model 

We have used MD calculations to study the influ- 
ence of surface roughness on liquid droplet contact angle 
and contact angle hysteresis. We have studied hydro- 
carbon liquid droplets on different self-affine fractal sur- 
faces. The nanodroplets contained 2364 octane molecules 
CsHig at T = 300 K, which is between the melting and 
boiling points of octane. The fractal surfaces were gen- 
erated by adding plane waves with random phases (see 
Sec. IID and Ref. [41(). Periodic boundary conditions are 
applied along the x and y directions. The periodically re- 
peated cell forms a rectangle L x x L y with L x = 506 A 
and L y = 75.9 A (see Fig- [9j) - The (non-contact) cylin- 
drical droplet diameter is about 104 A, and the size of the 
droplet-substrate contact area varies (for the hydropho- 
bic system) from w 115 A (case (a) in Fig. ITTJ]) to w 60 
A (case (c)). 

For most real surfaces usually there is a roll-off wave 
vector qo, below which the power spectrum C{q) of the 
surface is approximately constant. For q > q^ we assume 
the powe r sp ectrum has the power-law behavior C(q) ~ 
g -2(H+i) pjj corresponding to a self affine fractal surface 
with the fractal dimension D{ = 3 — H . Different fractal 
surfaces are obtained by changing the root mean square 
(rms) roughness amplitude cr, and the fractal dimension 
D{ . The roll-off wave- vector for the rough surface is (ft = 
27r/L with L = 38 A, and the magnitude of the short- 
distance cut-off wave vector q\ — n/a, where a = 2.53 
A is the substrate lattice constant. In the present work, 
we are mainly interested in how the rough structure of 
the substrate influence wetting. A curved upper wall has 
been used to speed up the formation of the droplet in the 



initial preparation [42j and to limit the gas-phase volume 
so that the droplet cannot (fully) evaporate. In some 
simulations we used a fiat upper wall in order to be able 
to push the droplet towards the substrate a bit more, 
in order to study the evolution of the receding contact 
angle. 

The lubricant molecules are described through the Op- 
timized Potential for Liquid Simulation (OPLS)[43l lil]: 
this potential is known to provide density and viscosity of 
hydrocarbons close to the experimental one. Each octane 
molecules comprises four units (particles), each particle 
corresponding to one chemical group CH 2 , CH 3 . The 
interaction between particles of different molecules is de- 
scribed by Lennard- Jones potentials. The intramolecular 
interactions include two body forces that keep the bond 
length C — C close to 1.53 A, three body forces imposing 
a preferred angle of 115° between the carbon atoms, and 
four body forces favoring a well defined torsion of the 
molecules. The four body forces apply to the sequence 
of carbon atoms C - C - C - C[45|. 

We used the Lennard-Jones interaction potential be- 
tween droplet atoms and substrate atoms. The L.-J. pa- 
rameters for a hydrophobic surface are chosen such that 
the Young contact angle is about 100° when a droplet 
sits on the flat surface. Because of the periodic boundary 
conditions and the size of our system, the liquid droplet 
forms a cylinder with the central line along the y-axis, 
as shown in Figures [5] and [TUJ We fit the density profile 
of the droplet to a cylinder (see Fig. ITTj) , and obtain the 
contact angle 6 = 103° for the droplet in contact with the 
flat hydrophobic substrate, while for the flat hydrophilic 
substrate 6 = 39 ±3°. 



B. Multiscale rough surfaces 

Many solid surfaces in nature, e.g., surfaces prepared 
by fracture (involving crack propagation), tend to be 
nearly self-affine fractal. Self-affine fractal surfaces have 
multiscale roughness, sometimes extending from the lat- 
eral size of the surface down to the atomic scale. A self- 
affine fractal surface has the property that if part of the 
surface is magnified, with a magnification which in gen- 
eral is appropriately different in the directory perpendic- 
ular to the surface as compared to the lateral directions, 
the surface "looks the same" [46j i.e., the statistical prop- 
erties of the surface are invariant under this scale trans- 
formation. 

The most important property of a randomly rough sur- 
face is the surface roughness power spectrum C(q). We 
assume that the statistical properties of the surface are 
translational invariant and isotropic so that C{q) depends 
only on the magnitude q = |q| of the wavevector q. For a 
self-affine surface the power spectrum has the power-law 
behavior C(q) ~ q~ 2i - H+1 \ where the Hurst exponent 
H is related to the fractal dimension Df = 3 — H. Of 
course, for real surfaces this relation only holds in some 
finite wave vector region qo < q < q±. Note that in 
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FIG. 10: Snapshots for different root mean square roughness, 
(a) the droplet is in contact with the flat substrate, (b) and (c) 
are for rough substrates with the root mean square amplitude 
a = 2.3 A and a = 4.8 A, respectively. Adapted from Ref. 

El. 
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FIG. 11: Determination of the contact ang le 8 for the flat 
substrate. Side view. Adapted from Ref. |ll| . 



many cases there is roll-off wavevector go below which 
C(q) is approximately constant. We have generated self 
affine fractal surfaces by adding plane waves with random 
phases and ap propria tely chosen weights, as described in 
detail in Ref. |4l|,l46| . 

In Fig. [12] we show the average slope £ and the ratio 
A/Ao between the surface area A and the nominal (or 
projected) Aq surface area, as a function of the root mean 
square roughness a when Hurst exponent H = 0.8, and 






FIG. 12: The average slope £ and the ratio A/Ao between the 
actual A and the nominal (or projected) Ao surface 3X6 EL j ELS 
a function of the root mean square roughness a when Hurst 
exponent H — 0.8, and as a function of Hurst exponent H for 
a = 3 A. 



FIG. 14: The contact angle 9 as a function of Hurst exponent 
H for the rms roughness a = 3 A. The circles and squares have 
the same meaning as that in Fig. 1131 The fractal dimension is 
D t = 3 - H. Adapted from Ref. [lj. 
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FIG. 13: The contact angle as a function of the root mean 
square roughness a. The circle points are numerical results 
from the simulations, while the square points are obtained 
from the Cassie model (see Eq. [7J). Each data point is an 
average over several snap-shot configurations. The fractal di- 
mension is Df = 2.2. Adapted from Ref. [Til ]. 
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FIG. 15: The contact angle 9 as a function of the normalized 
surface area A/Ao when the root-mean-square amplitude a 
increases for fixed Hurst exponent (H = 0.8, the solid circles) 
(from Fig. I13|l and when the Hurst exponent decreases for a 
fixed root-mean-square amplitude (cr = 3 A, the solid squares) 
(from Fig. Q3} . 



A. Static contact angle on hydrophobic surface 



as a function of Hurst exponent H for a = 3 A. 



IV. NUMERICAL RESULTS 

We present numerical results for the contact angle and 
contact angle hysteresis for both hydrophilic and hy- 
drophobic systems. The substrate surfaces are assumed 
to be self afHne fractal, but we have varied the fractal di- 
mension and the root-mean-square roughness amplitude. 



The apparent contact angle, 9q, as a function of the 
root mean square roughness (rms), is shown in Fig. [T3] 
with the fractal dimension Df = 2.2. There is a strong 
increase in #o with increasing rms-roughness amplitude. 
Fig. [JJ] shows how 6q depends on the Hurst exponent 
H = 3 — Df. Note that #o is almost independent of H. 

Accordingly to the Wenzel equation, the apparent con- 
tact angle 8 depends only on the surface roughness via 
the ratio r = A/Ao. Fig. fT2l shows that as H decreases 
from 1 to 0.4 (i.e., Df increases from 2 to 2.6), A/Aq in- 
creases by ~ 50%. However, the MD-calculations show 
that the apparent contact angle #o is almost independent 
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FIG. 16: (a) The height probability distribution for hy- 
drophobic surface both fiat (squares) and rough (circles), (b) 
The height probability distribution for hydrophilic surface 
both flat (squares) and rough (circles). 



of the fractal dimension, see Fig. [14] This is also illus- 
trated in Fig.Qjjjwhich shows the contact angle as a func- 
tion of the (normalized) surface area A/Aq for both cases. 
Thus the Wenzel equation cannot be used in the present 
situation. This is consistent with a visual inspection of 
the liquid-substrate interface which shows that on the 
rough substrates, the droplet is "riding" on the asperity 
tops of the substrate, i.e., the droplet is in a Cassie-likc 
state. In order to quantitatively verify this, we have cal- 
culated the distances u(x, y) between the bottom surface 
of the liquid drop and the rough substrate surface in the 
(apparent) contact area. From the distribution [47| 

P(u) = (S[u-u(x,y)}) 



of these distances [see Fig. [ToTa)] we obtain the fraction 
ip of the (projected) surface area where contact occurs: 



du P(u), 



where u\ is a cut-off distance to distinguish between con- 
tact and no-contact regions, which has to be comparable 
to the typical bond distance (we use u\ = 4 A). Note that 




FIG. 17: The advancing contact angle f? a evolution for hy- 
drophobic nanodroplet. f9 a is measured when the solid/liquid 
contact area increases. 



due to the thermal fluctuations tp — tpo for flat surface is 
less than 1. Using the normalized <f> = ip/ipo, the Cassie 
model (with r = 1) predicts the variation of the con- 
tact angle with a and H given in Fig. IT51 and [T4l (square 
points) . 

Fig. Q2] shows that the apparent contact angle 00 
increases strongly with increasing rms-roughness am- 
plitude, at fixed fractal dimension Df = 2.2, while 
it is nearly independent of the fractal dimension D[ 
(see Fig. [T4|) . Since increasing the fractal dimension 
at constant rms-roughness amplitude mainly increases 
the short- wavelength roughness, we conclude that the 
nanoscale wavelength roughness doesn't matter so much 
in determining the contact angle for hydrophobic sur- 
faces, while the long wavelength roughness plays an im- 
portant role. We attribute this fact to the strong ther- 
mal fluctuations in the height (or width) u of the liquid- 
solid interface which occur at the nanoscale even for the 
flat substrate surface. Note also that in our model the 
wall-wall interaction is long-ranged, decaying effectively 
as ~ 1/u 3 , so there will be a contribution to the in- 
terfacial energy also for non-contacting surfaces which, 
of course, is not rigorously included in the macroscopic 
Cassie model. 



Dynamic contact angle: 
hysteresis 



Contact angle 



The advancing contact angle 6* a is measured when 
the solid/liquid contact area increases, while the reced- 
ing contact angle 9 r is measured when the contact area 
shrinks. If the difference 9 a — 6 r is nonzero, the liquid- 
substrate system exhibits contact angle hysteresis. 




FIG. 18: The receding contact angle 8? evolution for hy- 
drophobic nanodroplet. 9 r is measured when the solid/liquid 
contact area shrinks. 
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FIG. 19: The advancing (circles) and receding (squares) con- 
tact angle 9, for hydrophobic substrate, as a function of time. 
The root-mean-square roughness of the substrate is rms = 4.8 
A. e = 4 meV and ro = 3.28 A. The thermal equilibrium con- 
tact angle has been reached after a few nanoseconds, irrespec- 
tive of whether the initial contact angle is larger or smaller 
than the equilibrium angle. 



1. Hydrophobic surfaces 

Figures [17] and [18] show the time-evolution of the ad- 
vancing contact angle 9 a and receding contact angle 9 r , 
respectively, for a nanodroplet on a rough, hydropho- 
bic substrate. The former has been obtained by placing 
the droplet close to the substrate, so that the drop will 
spontaneously spread under the adhesive interaction with 
the substrate. The contact angle evolves in time from 
9 = 180° (non-contact) in panel [TTT a) to its asymptotic 
value 9 a = 140°, reached after 3ns in panel [TWh) . 

The receding contact angle was simulated by squeezing 
the droplet into a pancake- like shape with the upper wall. 
The interaction between the atoms of the upper wall and 
drop atoms is given by the repulsive term of a Lennard- 
Jones potential, i.e. V(r) — 4eo(ro/r) 12 . The lack of at- 
traction with the top surface allows us to suddenly pull 
the wall up, leaving the drop in the configuration of panel 
ITST a). The free drop increases the contact angle up to 
the asymptotic value 9 r . Figure [T91 shows the time evo- 




Time (ns) 



150 



X 120 



60 



(b) o=3A 






J' 












7 




H=0.4 — 1 — 
H=0.5 — X— 
H=0.6 — X— 
H=0.7 — a— 
H=0.8 
H=0.9 

Efcl.O — •— 



1 

Time (ns) 



FIG. 20: The time evolution of the receding contact angle for 
the hydrophobic droplets, on the substrates with (a) various 
root-mean-square (rms) roughness while the Hurst exponent 
H = 0.8, and (b) various Hurst exponent H while the root- 
mean-square roughness (rms = 3 A). 



lution of the contact angle for these two cases. Note the 
strong time oscillations of the contact angle which are 
due to the rearrangement of the liquid molecules at the 
solid-liquid interface. The corresponding energy barri- 
ers are small compared to the thermal fluctuations (see 
Sec. IIC). However, after a few nanoseconds we find that 
both the receding and advancing contact angle fluctuate 
around the same average value; thus no contact angle 
hysteresis is observed for the hydrophobic system. 

To be sure that in any system there is no contact angle 
hysteresis, we performed extensive simulations on var- 
ious substrates with different root-mean-square rough- 
ness amplitudes (rms) (see Fig.l2"DTa)). and with different 
Hurst exponents H (see Fig.[20Tb)). The receding contact 
angle reaches its asymptotic value within about 2 or 3 
nanoseconds. In Figure [^01 one can see a relatively broad 
range of receding contact angles for substrates wih dif- 
ferent rms roughness. Conversely, substrates with differ- 
ent Hurst exponents show nearly the same contact angle. 
The root-mean-square roughness is mainly determined by 
long wavelength roughness of the surface. Increasing the 
fractal dimension Df signifies that the short wavelength 
roughness increases. Thus, one can see that the contact 
angle is more sensitive to the long wavelength roughness 
of the substrate than to the short wavelength roughness. 
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FIG. 21: Advancing contact angle evolution for hydrophilic 
nanodroplet. The root-mean-square roughness of the sub- 
strate is rms = 4.8 A. The energy parameter and the equilib- 
rium distance in L.-J. potential are e = 8 meV and ro = 3.28 

A. 



This agrees with the results in Fig. [13] and [14] 

A comparison of these results with those of the cor- 
responding simulations for the advancing contact angle, 
confirms that there is no hysteresis. This is in drastic 
contrast to simulation studies we have performed for hy- 
drophilic surfaces (see below), where surface roughness 
results in strong pinning of the boundary line; for such 
surfaces it is therefore impossible to study static droplet 
contact angles (as observed on macroscopic time scales) 
using molecular dynamics. 

Comparing the form of P(u) for the flat and the most 
rough surfaces shows that the system is in a Cassie-like 
state, but at the nanoscale the difference between the 
Cassie state and the Wenzel state is not so large due 
to the thermal fluctuations. That is, already for the flat 
surface strong thermal fluctuations at the interface result 
in nanosized regions where the separation between the 
solid and the liquid is much larger than the natural (low- 
temperature) binding separation. When the substrate is 
rough the fluctuations become even larger and the system 
is in a state which is more Cassie-like than Wenzcl-likc. 
This also explains why no hysteresis is observed: The 
Wenzel state is probably the (low temperature) energy 
minimum state (see Sec. IID), but squeezing the droplet 
into a pancake shape does not push the system perma- 
nently into the Wenzel state (where pinning effects may 
be very important) because even if it would go into this 
state temporarily, the free energy barrier separating the 
Cassie and Wenzel states is so small that thermal fluctu- 
ations would quickly kick it back to the Cassie-like state. 




FIG. 22: Receding contact angle evolution for hydrophilic 
nanodroplet. The root-mean-square roughness of the sub- 
strate is rms = 4.8 A. The energy parameter and equilibrium 
parameter in L.-J. potential are e = 8 meV and ro = 3.28 A. 
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FIG. 23: The advancing (circles) and receding (squares) con- 
tact angle 8, for the hydrophilic substrate, as a function of 
time. The root-mean-square roughness of the substrate is 
rms = 3 A, and the L.-J. substrate-liquid interaction param- 
eters e = 8 meV and ro = 3.28 A. 



2. Hydrophilic surfaces 

Let us now consider the case where the liquid droplet 
contact angle on the flat surface 9 < 90° (hydrophilic 
system). We choose the energy parameter and the equi- 
librium distance in L.-J. potential, associated with the 
liquid-solid atom interaction, as e = 8 meV, ro = 3.28 A 
respectively. This gives 6 w 70°. Figures EB [22] and [23] 
for hydrophilic droplet, are analogous to Figures [T71 [T81 
and [TUl respectively, for hydrophobic droplet. In Fig. |2"T1 
we show the time dependence of the advancing contact 
angle for the hydrophilic nanodroplet. The root-mean- 
square roughness of the substrate is rms = 4.8 A. 

In Fig. [HI we show the time evolution of the receding 
contact angle for the hydrophilic nanodroplet. 

Fig. [221 shows the advancing (circles) and receding 
(squares) contact angle 6 as a function of time, on a 
hydrophilic substrate with root-mean-square roughness 
rms = 3 A. Note that the thermal equilibrium contact 
angle cannot be reached on the time-scale of the sim- 
ulations. Note also that the fluctuation in the contact 




FIG. 24: The natural separation between the molecules in 
a liquid is usually considerably larger than the atom-atom 
separation on the substrate surface. This implies that the 
fluid molecules cannot "follow" the atomic scale roughness so 
that the fluid molecules will naturally be in a Cassie-like state 
with respect to the shortest substrate roughness components 
determined by the substrate nearest neighbor atom-atom sep- 
aration. 



angle are much smaller than for the hydrophobic system 
(Fig. [TH]). This is, of course, due to the fact that for 
the hydrophilic system the liquid-substrate interaction is 
much stronger, and the barriers for the rearrangement 
of liquid molecules at the substrate-liquid interface much 
higher than for the hydrophobic system. 

Finally, in Fig. [TWb'l we show the height probability 
distribution for the hydrophilic surface for both the flat 
(squares) and rough (circles) hydrophilic surface. Note 
that the fluctuations in the liquid-solid separation at the 
interface is much smaller on the hydrophilic surface than 
on the hydrophobic surface [ISf a)] . 



V. DISCUSSION 

In most practical cases it is not possible to modify the 
surface roughness without simultaneously affecting the 
chemical nature of the surface. While this is obvious for 
crystalline materials, where surface roughening will re- 
sult in the exposure of new lattice planes with different 
intrinsic surface energy, it may also hold for amorphous- 
like materials, where surface roughening may result in 
a more open atomic surface structure, with an increased 
fraction of (weak) unsaturated bonds. In our model study 
a similar effect occurs, and some fraction of the change in 
the contact angle with increasing root-mean-square am- 
plitude may be associated with this effect. However, the 
most important result of our study, namely that the con- 
tact angle is mainly determined by the long-wavelength 
roughness, should not be affected by this fact. 

Another reason for why the short-wavelength (atomic) 
roughness may influence the liquid contact state differ- 
ently from long-wavelength roughness has to do with the 
fact that the natural separation between the molecules 
in a liquid is usually considerably larger than the atom- 
atom separation on the substrate surface. This implies 
that the fluid molecules cannot "follow" the atomic scale 
roughness (see Fig. so that the fluid molecules will 



naturally be in a Cassie-like state with respect to the 
shortest substrate roughness wavelength components, de- 
termined by the substrate nearest neighbor atom-atom 
separation. 

It is important to note that the discussion in this pa- 
per is also relevant for the contact between macroscopic 
liquid drops and rough substrates. That is, if the solid- 
liquid interface is studied with nanometer resolution A, 
then the strong fluctuations (in time and space) at the 
liquid-solid interface discussed above will also be ob- 
served for a macroscopic droplet, and the interfacial en- 
ergy 7si(C) for C = D/X (where D is the diameter of the 
nominal contact area between the (macroscopic) liquid 
droplet and the substrate) will be the same as obtained 
above for nanodroplets at the same resolution A. 



VI. SUMMARY AND CONCLUSION 

We have discussed under which condition the Wenzel 
or Cassie state is favorable on randomly rough surfaces. 
We performed molecular dynamics simulations to study 
contact angle and the contact angle hysteresis on hy- 
drophobic and hydrophilic surfaces. The contact angle 
on hydrophobic surfaces depends strongly on the root- 
mean-square roughness of the substrate, but is nearly 
independent of the fractal dimension. For hydrophobic 
surfaces, there is no contact angle hysteresis due to strong 
thermal fluctuations at the nanoscale. For hydrophilic 
surfaces we observe contact angle hysteresis due to pin- 
ning effects resulting from the much higher energy bar- 
riers for rearrangement of liquid molecules at the solid- 
liquid interface. This indicates that on randomly rough 
hydrophobic surfaces the Cassie-like state often prevails, 
at least for nanoscale droplets. We have found that ther- 
mal fluctuations play an important role at the nanoscale, 
which leads to the enhanced hydrophobicity by surface 
roughness. It is of particular importance to design and 
build super hydrophobic surfaces. 



VII. APPENDIX: DISTRIBUTION OF 
SURFACE SLOPES FOR RANDOMLY ROUGH 
SURFACES 

In this appendix we derive the distribution of sur- 
face slopes for randomly rough surfaces. We discuss un- 
der which conditions one expects the Wenzel and Cassie 
states to prevail. 

The surface area A and the average surface 
slope £ 

Consider a randomly rough surface and let h(x) de- 
note the height profile measured from the average plane 
so that (h(x)) = 0, where (..) stands for ensemble averag- 
ing, or (equivalently) averaging over the surface area. We 
assume that /i(x) is a Gaussian random variable charac- 



terized by the power spectrum 
1 



C(q) = 



(2tt)' 

Note that if we write 



d 2 x (h(x)h(p))t 



- / q x 



h(x) = J d 2 q /i(q)e lq - : 



where 



then 



h(q) 



(2nf 



d 2 x /i(x)e- ,q - > 



(h(q)h(q')) = S(q + q')C(q). 
Sometimes it is also convenient to use 



(/i(q)M-q)) = 



An 



(2n) 



;C(q), 



(Al) 



(A2) 



where Aq is the surface area. In deriving (A2) we have 
used that 



*(q-q) 



1 



(2tt)2 



d 2 x e i( *-ti-= = 



An 



(27T) 



If the surface roughness amplitudes h(q) are assumed 
to be Gaussian random variables, one can show that the 
(normalized) surface area[Hj 



.4 
T 



, =^-= I dx ii+xef' 2 ' 



where 



d 2 q q 2 C(q) = 2ir dq q 3 C(q) 



Let us calculate the rms surface slope. We get 
((Vh) 2 ) = / d 2 qd 2 q' (iq a )(iq' a )(h(q)h(q'))< 



d 2 q q 2 C(q) = £ 2 



«(q+q') 



Using (Al) this gives 
((S7h) 2 ) 



Thus, for a Gaussian random surface both the average 
slope and the increase in the surface area is determined 
by the parameter £. For non-random surfaces this is no 
longer the case. 

Surface slope probability distribution 

Let h(x, C) denote the height profile after having 
smoothed out surface roughness with wavelength shorter 
than A = L/C,. For example, we may define 

fc(x,C)= / d 2 qh(q)e^, 
J q<qi 



where q\ = qi,C (where q^ — 2tt/L). We will refer to £ as 
the magnification. Thus, when we study the surface at 
the magnification £ we will only detect surface roughness 
with wavelength components larger than A = L/C,. 

We will now derive an equation of motion for the sur- 
face slope probability distribution function 

P( S ,C) = («5( S -V/i(x,0)). 

We assume that the surface roughness amplitudes h(q) 
are independent random variables. In this case, if we 
write 



we get 



/i(x, C + SC.) = h(x, C) + Sh, 



But 



= / d 2 s'{5{s' - VSh)){S{s - s' - Vh{x, ())) 



J d 2 s'(S(s' - VSh))P(s - s', C). (A3) 



(5(s< - V5h)) = -J-p J d 2 q(e^'~™) ) 



^Jd 2 q (l-i((q.V^) 2 ) 



2\ \ giq-s' 



(2 



= 8(s') + \((V a Sh) (V^ffi)^^(s'). 

Substituting this result into (A3) and expanding the left- 
hand-side to linear order in SC gives 

f(^) = ^((v^)(v^))A_^p (s , c) . m 



But 



((V a Sh) (VpSh)} = / d 2 qd 2 q' 
JqU 



x{k a )(iq'f))(h(q)h{q')) 



M\ e (q+q')-x 



9i(C+-5C) 



d 2 q q a q/3C(q) 



1 f qL(C+SC) 

-A, 1 



d'q q'C(q) 



qU 



Thus 



ir5 af 3q L 6(q 3 C(q). 



Substituting this result in (A4) gives the following 
diffusion- like equation for P(s, £): 

BP 



(A5) 



where 



_d d_ 

ds a 9s a 



where the "diffusivity" 



D(C) = ^q L q 3 C(q), 



(A6) 



with q = qi,C' 



Solution of the diffusion equation 

The function P(s, £) describes the probability to ob- 
serve a surface slope or gradient s = V/i(x) when the sys- 
tem is studied at the magnification £. When the system 
is studied at the lowest magnification £ = 1 the surface 
appears flat and smooth so that the gradient vanishes, 
i.e., 

P( S ,l) = <<5(s-VMx,l))>=£( S ). (AT) 

We also require that there is no infinite high slopes, i.e., 

P(s,C)^0 as |s|-+oc. (A8) 

Let us determine the solution to (A5) which obeys the 
"initial" condition (A7) and the boundary condition 
(A8). ft is clear that the solution is given by 



P(a,C) = -L e 

7TSf 



(A9) 



where the width parameter s\(Q depends on the magni- 
fication: 



4 I dCD(C) 



2tt 



dqq 3 C(q)=e(0- 



Note that P is normalized, 



d 2 s P(s,C) = l, 



and that the width of the Gaussian distribution P in- 
creases with increasing resolution, i.e., when the surface 
is studied at higher and higher resolution, steeper and 
steeper surface slopes will be detected. 



Distribution function P(sq,() 




In what follows we will need the fraction P(s ,C) 01 
the total surface area where the slope s < sq, where Sq 
is a fixed number between zero and infinite. Let A±(() 
be the surface area, projected on the xy-plane, where the 
surface slope s < s . We then have 



P(s ,O=A±(O/A . 



Using the definition 



p(s,c) = (5(s-vMx,C))} 



j- J d 2 x a(s-v/i(x,o) 



we get 



P(s ,C)= / d 2 sP(s,C). 

J s<s 

Using (A9) this gives 

P( So ,C) = i-^ (so/si)2 . 



(A10) 



Surface area with slope below tan0 

Consider a liquid in contact with a rough substrate. 
The contact angle 9 is determined by Young's equation: 



cos# = 



7sv - 7si 
7iv 



where j a i, 7 SV and 7i v are the solid- liquid, solid- vapor and 
liquid- vapor interfacial energies, respectively. Note that 
if n is the normal to the solid surface and z the normal to 
the liquid surface, which we assume to be parallel to the 
average surface plane (see Fig. [55]), then cos9 = — z n. 
Since 



n = 



(-VM) 

(I + (V/j) 2 ) 1/2 



we get 



costf = - (I + (V/i) 2 ) 1/2 



|tan0| = \Vh\. 



Thus, using (AIO) the fraction of the surface where the 
surface slope is below tan6> is 



P(tan0) = I 



-(tanfl/£) 2 



(All) 



FIG. 25: A fluid in contact with a rough substrate. The 
contact angle 9 is determined by Young's equation. 
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